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A  comprehensive  data  set  from  the  ocean  and  atmosphere  was 
obtained  just  north  of  the  Monterey  Bay  as  part  o I  the  Mon¬ 
terey  Bay  2006  (MB06)  field  experiment.  The  wind  stress,  heat 
fluxes,  and  sea  surface  temperature  were  sampled  by  the  Naval 
Postgraduate  School's  TWIN  OTTER  research  aircraft.  In  situ  data 
were  collected  using  ships,  moorings,  gliders  and  AUVs.  Four 
data-assimilating  numerical  models  were  additionally  run,  includ¬ 
ing  the  Coupled  Ocean/Atmosphere  Mesoscale  Prediction  System 
(COAMPS®)  model  for  the  atmosphere  and  the  Harvard  Ocean 
Prediction  System  (HOPS),  the  Regional  Ocean  Modeling  System 
(ROMS),  and  the  Navy  Coastal  Ocean  Model  (NCOM)  for  the  ocean. 

The  scientific  focus  of  the  Adaptive  Sampling  and  Prediction 
Experiment  (ASAP)  was  on  the  upwelling/relaxation  cycle  and 
the  resulting  three-dimensional  coastal  circulation  near  a  coastal 
promontory,  in  this  case  Point  Afto  Nuevo,  CA.  The  emphasis  of  this 
study  is  on  the  circulation  over  the  continental  shelf  as  estimated 
from  the  wind  forcing,  two  ADCP  moorings,  and  model  outputs.  The 
wind  stress  during  August  2006  consisted  of  3-10  day  upwelling 
favorable  events  separated  by  brief  1  -3  day  relaxations  During  the 
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first  two  weeks  there  was  some  correlation  between  local  winds 
and  currents  and  the  three  models’  capability  to  reproduce  the 
events.  During  the  last  two  weeks,  largely  equatorward  surface 
wind  stress  forced  the  sea  surface  and  barotropic  poleward  flow 
occurred  over  the  shelf,  reducing  model  skill  at  predicting  the  cir¬ 
culation.  The  poleward  flow  was  apparently  remotely  forced  by 
mesoscale  eddies  and  alongshore  pressure  gradients,  which  were 
not  well  simulated  by  the  models.  The  small,  high -resolution  model 
domains  were  highly  reliant  on  correct  open  boundary  conditions 
to  drive  these  larger-scale  poleward  flows.  Multiply-nested  mod¬ 
els  were  no  more  effective  than  well-initialized  local  models  in  this 
respect. 

€  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  Adaptive  Sampling  and  Prediction  Experiment  (ASAP)  was  launched  during  the  summer  of 
2006  with  the  overall  goal  to  learn  how  to  deploy,  direct,  and  utilize  autonomous  vehicles  most  effi¬ 
ciently  to  sample  the  ocean,  assimilate  the  data  into  numerical  models  in  real  or  near-real  time,  and 
predict  future  conditions  with  minimal  error.  A  fleet  of  four  Spray  (Davis  et  al.,  2002)  and  six  Slocum 
gliders  under  autonomous  control  were  deployed  in  late  July  to  cooperatively  sample  a  control  volume 
in  an  upwelling  center  for  a  period  of  30  days.  The  scientific  focus  was  to  close  the  heat  budget  in  a 
three-dimensional  upwelling  center  and  thus  illuminate  the  along-  and  across-shore  transport  and 
fundamental  dynamics  of  the  process.  The  glider  fleet  patrolled  the  boundaries  and  interior  of  the  con¬ 
trol  volume  and  a  research  aircraft  observed  the  fluxes  through  the  sea  surface.  Two  bottom-mounted 
acoustic  Doppler  current  profilers  (ADCPs)  were  deployed  on  the  54  and  92  m  isobaths  to  monitor 
the  currents  over  the  continental  shelf  and  provide  ground  truth  for  the  newer,  more  innovative 
instrumentation. 

Based  on  previous  work,  the  evolution  of  coastal  upwelling  near  Point  Ano  Nuevo  can  be  described 
as  follows:  By  comparing  aircraft-sensed  winds  and  SST,  Ramp  et  al.  (2009)  determined  that  cold 
water  <13°C  first  appears  at  the  Point  within  12-24h  of  the  onset  of  upwelling-favorable  winds.  As 
the  wind  continues  to  blow  for  5-10  days,  a  cold  plume  forms  which  may  move  further  offshore  or  turn 
southward  across  the  Monterey  Submarine  Canyon  (MSC)  (Rosenfeld  et  al.,  1994:  Ramp  et  al.,  2005; 
Haley  et  al.,  2009).  A  warm,  fresh,  anticyclonic  feature  (the  “Monterey  Bay  Eddy”  or  MBE)  often  exists 
offshore,  which  interacts  with  the  cold  plume  in  some  complex  way  (Rosenfeld  et  al.,  1 994;  Paduan  and 
Rosenfeld,  1 996;  Ramp  et  al.,  2005).  It  has  now  been  determined  during  the  second  Autonomous  Ocean 
Sensing  Network  (AOSN-U)  experiment  that  the  position  and  strength  of  the  anticyclonic  feature  can 
determine  the  bifurcation  of  the  upwelling  front  (Ramp  et  al.,  2009).  When  the  MBE  lies  offshore,  the 
upwelling  plume  is  transported  offshore  by  Ekman  advection.  When  the  feature  has  moved  inshore, 
the  combination  of  blocking  and  geostrophic  current  force  the  plume  to  move  southward  across  the 
mouth  of  the  bay.  Time  series  of  sea  surface  temperature  from  aircraft  and  AVHRR  images  show  that 
the  MBE  moves  rapidly  onshore  during  wind  relaxation  events  and  retreats  back  offshore  when  the 
winds  re-intensify  (Rosenfeld  et  al.,  1994;  Ramp  et  al.,  2005).  Moored  observations  however  have 
also  shown  an  onshore  and  southward  translation  of  the  MBE  without  a  wind  relaxation  (Ramp  et  al., 
1997a,b).  The  mechanism  for  this  across-shore  eddy  translation  is  thus  not  fully  understood.  How  far 
the  eddy  forcing  penetrates  up  onto  the  continental  shelf  has  also  not  been  well  observed.  Since  these 
features  are  more  than  1000  m  deep  (Ramp  etal.,  1997a,b)one  expects  that  due  to  rotational  rigidity, 
their  penetration  onto  the  shelf  would  be  limited,  similar  to  Gulf  Stream  Warm  Core  Rings  (Ramp, 
1989). 

Much  less  obvious  is  the  ocean  dynamics  during  the  wind  relaxation  events,  defined  here  to  be 
periods  of  poleward  or  no  wind.  During  these  events,  barotropic  poleward  flows  typically  appear  over 
the  continental  shelf  along  the  west  coast  of  the  U.S.  Observational  examples  of  this  include  the  Coastal 
Ocean  Dynamics  Experiment  (CODE)  region  off  Point  Arena,  CA  (Winant  et  al.,  1987;  Lentz,  1987), 
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further  south  off  Point  Sur,  CA  (Chelton,  1984;  Ramp  and  Abbott,  1998),  and  south  of  Cape  Blanco, 
OR  (Ramp  and  Bahr,  2006).  Recent  modeling  studies  attribute  this  current  to  a  poleward  alongshore 
pressure  gradient  force  (PGF)  south  of  the  major  headlands,  which  encourages  both  poleward  flow 
near  the  coast  and  onshore  transport  in  geostrophic  balance  with  the  poleward  PGF  (Gan  and  Allen, 
2002a, b).  Since  these  gradients  are  typically  remotely  forced  over  very  large  scales,  there  is  usually  less 
coherence  between  local  winds  and  currents  during  relaxation  events  than  during  upwelling  events, 
making  them  much  more  difficult  to  reproduce  with  local,  high-resolution  numerical  models.  In  fact, 
uncertainties  of  coherent  structures  estimated  by  model  simulations  (Lermusiaux,  2006;  Lermusiaux 
et  ah,  2006)  during  such  relaxation  periods  are  large. 

In  this  paper,  a  rigorous  description  of  the  currents  over  the  continental  shelf  just  south  of  Point 
Ano  Nuevo  and  their  response  to  local  and  remote  forcing  is  presented  using  in  sifu  data  from  moored 
ADCPs,  remotely  sensed  data  from  a  low-flying  aircraft,  and  model  output  from  three  numerical  models 
run  during  the  ASAP  experiment.  The  paper  is  intended  to  set  the  stage  for  a  more  detailed  dynamical 
analysis  of  the  three-dimensional  circulation  in  the  ASAP  region  using  all  available  assets.  Follow¬ 
ing  a  basic  description  of  the  observational  tools  and  techniques  in  Section  2,  the  observed  winds 
and  currents  are  described  in  Section  3  and  compared  using  correlation  and  spectral  analysis.  The 
three  models,  including  their  grids,  initialization  and  data  assimilation  techniques,  and  computational 
schemes  are  set  forth  in  Section  4,  where  the  model  outputs  are  also  quantitatively  compared  with  the 
observed  currents  at  the  two  ADCP  locations.  The  results  are  discussed  in  Section  5  and  summarized 
in  Section  6. 

2.  Data  and  methods 

2./.  Real-time  ADCP  data 

The  moorings  deployed  for  the  August  2006  experiment  consisted  of  two  300  kHz  acousti  c  Doppler 
current  profilers  (ADCPs)  in  trawl-resistant  bottom  mounts  (TRBMs)  located  6.55  km  apart.  ADCP 
1  was  located  at  36C55.336'N,  122  07.344  W,  54  m  depth,  and  ADCP  2  was  located  at  36  53.683'N, 
1 22°  1 1 ,244'W,  92  m  depth  ( Fig.  1 ).  The  data  were  sampled  in  4-m  bins  from  the  bottom  to  the  surface, 
however  the  uppermost  3  bins  ( 1 2  m)  were  lost  due  to  side-lobe  reflection.  Bottom-mounted  moorings 


1 22.5  122.4  122.3  122.2  122.1  122.0  121.9  121.8 


Fig.  1.  Map  of  the  greater  Monterey  Bay  region  showing  the  locations  of  the  ADCPs  and  wind  buoys.  On  the  continental  shelf, 
the  20,  50,  and  1 00  m  contours  are  shown. 


S.R.  Ramp  et  ai  /  Dynamics  of  Atmospheres  and  Oceans  52  (201 1 )  192-223 


195 


are  advantageous  because  they  have  no  extraneous  movement,  and  allow  faster  sampling  and  less 
averaging  to  achieve  a  good  signal-to-noise  ratio  than  a  surface  mooring.  The  basic  1-min  time  series 
were  obtained  by  ensemble  averaging  60  1-s  pings  in  the  instrument.  Subsequent  averaging  and 
filtering  to  separate  the  frequency  bands  of  interest  were  accomplished  in  the  laboratory.  An  inverse- 
Fourier  truncation  method  was  used  to  remove  internal  waves  with  periods  of  5  h  or  less.  The  residual 
time  series  were  filtered  again  with  the  half-power  point  at  33  h  (Beardsley  et  al.,  1985)  to  separate 
the  tidal  and  inertial  motions  from  the  wind  forced,  mesoscale,  and  mean  currents.  A  principal  axis 
technique  was  used  to  rotate  the  u  and  v  components  into  across-  and  alongshelf  components.  This 
method  of  rotation  ensures  that  a  minimum  of  the  alongshelf  flow  will  be  translated  into  the  across- 
shelf  flow  (Kosro,  1 987).  The  principal  axis  rotation  agreed  with  the  local  orientation  of  the  bathymetry 
and  yielded  angles  of  -58  for  ADCP  1  and  -48n  for  ADCP  2,  with  respect  to  true  north. 

A  novel  aspect  of  the  experiment  was  that  the  data  were  transmitted  back  to  the  laboratory  in 
real-time  using  Benthos  Telesonar  acoustic  modems  operating  on  a  Seaweb  undersea  network  (Rice 
and  Green,  2008).  This  system  allowed  the  ADCPs  to  be  bottom-mounted  in  TRBMs,  completely  free 
of  surface  motion,  but  still  able  to  transmit  the  data  to  the  laboratory  in  real  time.  The  ADCPs  were 
programmed  to  send  back  the  most  recent  estimate  via  satellite  link,  every  1 0  min.  A  redundant  system 
with  two  surface  gateway  buoys  and  three  intermediate  relay  points  was  used  to  ensure  reliability  and 
keep  the  acoustic  power  output  levels  low.  The  data  were  moved  from  the  gateway  buoys  to  the  Naval 
Postgraduate  School  via  the  iridium  satellite  communications  network.  The  system  worked  flawlessly 
for  1 8  days  but  then  failed  due  to  (a)  sea  water  incursions  into  the  iridium  antennas  and  (b)  low  power 
on  the  surface  buoys.  The  real-time  current  data  were  used  to  understand  and  adapt  the  sampling 
paths  of  the  ASAP  gliders  and  AUVs. 


22.  Wind  observations 

The  wind  time  series  was  recorded  by  buoy  M2  located  at  36  42'N,  122D23'W  at  the  mouth  of 
the  Monterey  Bay  and  maintained  by  the  Monterey  Bay  Aquarium  Research  Institute  (MBAR1).  Buoy 
M2  was  slightly  farther  away  from  the  two  ADCPs  than  similar  buoy  Ml,  but  was  less  sheltered  by 
coastal  topography  and  therefore  more  representative  of  the  offshore  winds.  The  wind  time  series 
were  low-pass  filtered  with  the  same  33-h  filter  as  the  ADCP  data  and  rotated  22°  into  the  along-  and 
across-shore  components,  aligned  with  the  coastal  mountain  range.  This  is  a  different  rotation  than 
was  used  for  the  currents  to  match  the  local  bottom  topography. 


2.3.  Aircraft  data 

The  airborne  measurements  during  ASAP  were  obtained  using  the  TWIN  OTTER  aircraft  owned 
and  operated  by  the  Center  for  Interdisciplinary  Remotely-Piloted  Aircraft  Studies  (C1RPAS)  at  the 
Naval  Postgraduate  School.  Thirteen  flights  were  conducted  daily  between  August  1  and  15,  2006 
(except  August  6  and  13)  at  an  altitude  of  nominally  33  m  above  the  sea  surface.  The  low  alti¬ 
tude  was  chosen  to  pass  beneath  the  stratus  cloud  deck  that  often  covers  the  Monterey  Bay  in 
summer,  and  to  minimize  the  vertical  divergence  in  the  heat  and  momentum  flux  observations. 
The  flight  path  covered  the  area  from  the  center  of  the  bay  north  to  about  Pigeon  Point,  out  to 
35  km  offshore.  This  took  about  2.5  h  to  execute  beginning  usually  about  10  am  local  time,  subject 
to  weather  conditions.  To  collect  air/sea  fluxes  at  high  temporal  and  spatial  resolution,  the  sam¬ 
ple  area  focused  on  the  region  surrounding  Point  Ano  Nuevo  where  the  upwelling  plume  was  first 
developing.  The  atmospheric  parameters  sampled  by  the  aircraft  included  air  temperature  and  dew 
point,  atmospheric  pressure,  wind  speed  and  direction,  the  turbulent  fluxes  of  heat  and  momen¬ 
tum,  and  total  particle  number,  as  well  as  aerosol  and  cloud/fog  droplet  size  distributions.  In  the 
ocean,  sea  surface  temperature  (SST),  ocean-leaving  radiance  at  193  wavelengths,  sun  glint,  and  sur¬ 
face  roughness  were  observed.  The  SST,  air  temperature  and  dew  point,  and  wind  data  were  made 
available  for  model  assimilation  in  near-real  time,  while  the  other  parameters  required  additional 
post-processing. 
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3.  Observational  results 

3.1.  Atmospheric  forcing 

August  2006  was  a  fairly  typical  summer  month  off  central  California  with  periods  of  sustained 
upwelling  favorable  winds  at  about  10ms  1  separated  by  brief  periods  of  weak  poleward  or  no  winds 
(Fig.  2).  At  buoy  M2,  upwelling-favorable  wind  events  with  approximately  equal  strength  took  place 
during  August  1-3,  8-12, 16-25,  and  28-31.  Poleward  winds  less  than  2.5  ms  1  blew  during  August 
4-7,  13-15,  and  26-27.  A  unique  strength  of  the  ASAP  data  set  was  the  spatial  distribution  of  the 
wind  stress  as  sampled  by  the  TWIN  OTTER  aircraft.  The  wind  stress  was  computed  from  the  bulk 
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Fig.  2.  Time  series  of  (a)  wind  stress  from  MBAR!  buoy  M2,  (b)  across-shore  current  at  ADCP  1,  (c)  alongshore  current  at  ADCP 
1 ,  (d)  across-shore  current  at  ADCP  2,  and  (e)  alongshore  current  at  ADCP  2.  The  scale  of  they-axis  is  the  same  for  both  ADCPs, 
with  the  difference  in  plot  height  reflecting  the  difference  in  water  depth.  The  velocity  scales  are  different  for  the  across-  and 
alongshore  current  components  to  better  illustrate  the  much  smaller  across-shore  flow.  A  positive  sign  indicates  poleward 
and  onshore  flow.  The  vertical  gray  bars  on  the  wind  stress  plot  indicate  the  time  of  the  TWIN  OTTER  over-flights  The  vectors 
indicate  the  direction  the  wind  is  blowing  towards. 
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aerodynamic  flux  formulae  using  the  wind  speed,  air  temperature,  relative  humidity,  air  pressure, 
sea  surface  temperature,  and  altitude  as  inputs  (Fairall  et  al.,  1996).  A  turbulent  boundary  layer  was 
assumed  and  Monin-Obukov  similarity  theory  was  used  to  adjust  the  stress  down  to  10  m  height. 
The  assumption  was  validated  by  the  bulk  Richardson  Number  (R,*)  which  was  small  and/or  negative 
throughout,  indicating  a  fully  turbulent  boundary  layer.  The  only  exception  occurred  on  August  5,  an 
exceptionally  calm  day,  but  R,-  was  still  less  than  10  and  the  same  method  was  used  to  compute  the 
stress. 

A  selection  of  8  of  the  15  flights,  4  during  a  relaxation  event  and  four  during  upwelling,  show 
significant  sub-mesoscale  variability  in  the  offshore  wind  stress  distributions  (Fig.  3).  The  first  wind 
relaxation  event  actually  lasted  longer  near  Point  Ano  Nuevo  than  indicated  by  the  time  series  at  Buoy 
M2,  and  spanned  6  days  from  2  to  7  August.  Especially  noteworthy  during  this  time  was  the  wind  stress 
maximum  located  just  to  the  north  of  the  Point.  This  feature  looks  like  an  atmospheric  expansion  fan 
similar  to  those  previously  observed  in  the  lee  (to  the  south)  of  major  capes  during  upwelling-favorable 
winds  (Dorman  et  al.,  1995, 2000;  Gan  and  Allen,  2005;  Gan  et  al.,  2005),  however  such  features  have 
not  been  directly  observed  during  poleward  wind  events.  The  maximum  wind  stresses  in  these  features 
(0.08  N  m  2)  on  August  3rd  and  4th  were  actually  much  stronger  than  indicated  by  the  offshore  buoys 
and  were  counter  to  the  generally  held  idea  that  the  winds  are  strongest  offshore  and  weaken  as  the 
coast  is  approached.  The  weakest  winds  during  the  relaxations  were  near  the  mouth  of  the  Monterey 
Bay. 

During  the  upwelling  event  starting  on  August  8th,  all  aircraft  observations  showed  an  alongshelf 
wind  stress  maximum  off  Point  Ano  Nuevo  (Fig.  3).  This  maximum  was  strongest  on  August  8th  and 
11th  (>0.16  Nm  2)  and  slightly  weaker  on  the  9th  and  10th.  This  clear  maximum  in  the  alongshore 
wind  stress  may  contribute  to  the  formation  of  the  upwelling  center  off  Ano  Nuevo.  Consistent  with 
previous  studies,  there  was  a  wind  shadow  behind  the  coastal  mountains  in  the  northeast  corner  of 
the  Monterey  Bay  (Graham  and  Largier,  1997;  Ramp  et  al.,  2005).  This  lead  to  a  very  strong  gradient 
in  the  wind  stress  extending  southwest  from  Santa  Cruz. 

To  quantify  these  gradients,  the  vertical  component  of  the  wind  stress  curl  k  (Vxt)  was  computed 
on  the  same  grid  as  the  wind  stress  itself  (Fig.  4).  A  planer  least-squares  fit  was  used  over  a  circle  with 
radius  15  km  to  estimate  the  curl.  During  the  weak  poleward  winds  on  August  4-5,  the  curl  changed 
sign  down  the  middle  of  the  observing  region  about  20  km  offshore  and  was  positive  (upwelling 
favorable)  offshore  and  negative  (downwelling  favorable)  nearshore.  This  may  have  contributed  to 
the  rapid  warming  of  the  sea  surface  (>3°C)  between  the  4th  and  the  5th  off  Ano  Nuevo.  On  August 
8, 10,  and  11,  there  was  a  region  of  high  positive  curl  extending  southward  across  the  Monterey  Bay 
from  Santa  Cruz,  along  the  boundary  between  the  unobstructed  offshore  wind  and  the  wind  shadow 
in  the  inner  bay.  Using  the  conventional  formula  for  Ekman  pumping 


we 


k-(V  xt) 

Pf 


the  upwelling  velocity  we  was  computed  where  p  is  the  water  density, /is  the  Coriolis  parameter, 
and  the  observed  curl  maximum  =  6  x  10  6Nm~3.  The  result  =  5.8  md'1  and  was  comparable  to  the 
upwelling  due  to  coastal  divergence  at  the  same  time,  which  was  in  the  range  of  5-10 md  1.  This 
indicates  that  the  cold  water  extending  southward  across  the  Monterey  Bay  during  upwelling  events 
may  result  from  the  local  wind  stress  curl  maximum  in  addition  to  the  southward  advection  of  cold 
water  from  the  north.  This  is  consistent  with  previous  results  on  larger  spatial  scales  that  also  show 
that  the  vertical  velocities  due  to  Ekman  pumping  are  comparable  to  those  from  Ekman  transport  near 
headlands  along  the  central  California  coast  (Pickett  and  Paduan,  2003). 


32.  ADCP  currents 


The  most  obvious  feature  in  the  observed  currents  was  the  preponderance  of  poleward  flow  at  both 
locations  (Fig.  2).  One  particularly  strong  poleward  event  started  on  August  3rd  at  ADCP  2  and  August 
4  at  ADCP  1.  The  event  was  nearly  barotropic  at  both  moorings  with  a  slight  mid-depth  maximum 
(near  40cm s  and  ended  around  August  7.  The  temporal  agreement  with  the  August  4-7  wind 
relaxation  appears  coincidental,  in  that  poleward  flows  with  this  magnitude  did  not  appear  during 
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Rg.3.  Time  series  of  the  aircraft-observed  surface  wind  stress  during  (a-d)  the  August  2-5  wind  relaxation  event  and  (e-h)  the  August  8-12  upwelling  favorable  event.  The  white  dots 
indicate  the  ASAP  mooring  locations  (nearshore)  and  MBAR1  buoy  M2  (offshore).  The  units  on  the  color  bar  are  in  N  m  2.  ( For  interpretation  of  the  references  to  color  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  the  article.) 
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the  subsequent  wind  relaxation  events.  Plan-view  displays  of  the  model  output,  discussed  m  detail  in 
a  later  section,  gives  some  indication  that  this  strong  poleward  event  was  confined  to  the  continental 
shelf  and  was  associated  with  a  pulse  of  warm  water  exiting  the  Monterey  Bay.  There  was  likewise 
evidence  in  the  models  that  the  equatorward  flow  at  ADCP  2  centered  on  August  11-12  was  due  to 
this  instrument  being  located  in  the  coastal  upwelling  jet,  while  ADCP  1  was  located  well  shoreward 
of  the  jet  where  the  flow  remained  poleward. 

Two  subsequent,  weaker  poleward  flow  events  took  place  from  August  19-21  and  August  26-28 
with  no  obvious  phase  lag  between  moorings.  These  events  were  barotropic  in  nature  but  had  stronger 
currents  near  the  bottom  (order  1 5  cm  s'1 )  than  the  surface  (order  5- 1 0  cm  s*1  )•  These  poleward  flows 
took  place  during  periods  of  no  winds  or  weak  equatorward  winds  and  thus  appear  to  be  wind-related. 
The  vertical  current  profiles  during  these  times  could  be  modeled  as  the  sum  of  a  barotropic  poleward 
pressure  gradient  force  over  the  entire  water  column  and  weak  equatorward  Ekman  forcing  in  the 
upper  20  m.  This  Ekman  flow  was  clearly  observed  in  the  across-shore  component  as  offshore  flow 
(blue  colors  in  Fig.  2)  during  the  stronger  upwelling-favorable  wind  events,  as  for  instance  during 
August  10-11,  16-17,  22-25,  and  29-31,  most  clearly  at  ADCP  2.  In  general,  ADCP  2  showed  more 
visual  correlation  with  the  surface  wind  stress  than  ADCP  1,  and  this  was  supported  by  the  spectral 
analysis  (not  shown).  Significant  spectral  peaks  were  observed  for  periods  of  10  days,  5  days,  and  the 
tides.  A  coherence  and  phase  analysis  showed  that,  on  average,  the  5-1  Od  winds  lead  the  alongshore 
currents  by  1 7.9  h  at  both  instruments.  This  is  close  to  the  local  inertial  period  (20  h)  and  represents  the 
spin-up  time  for  the  alongshore  flow  under  wind  forcing.  At  ADCP  2,  the  across-shore  currents  at  the 
near-surface  and  near-bottom  were  coherent  with  the  alongshore  wind  stress  with  opposite  phase 
( 1 51  °).  consistent  with  Ekman  dynamics.  At  ADCP  1,  there  was  no  coherence  between  alongshore 
winds  and  across-shore  currents  at  any  depth.  This  may  be  due  to  the  coastal  constraint,  as  Ekman 
turning  is  reduced  the  closer  the  location  is  to  the  coast  (Ekman.  1905;  Ramp  and  Abbott.  1998). 

Also  obvious  in  the  plots  are  the  much  shorter  space  and  time  scales  in  the  across-shore  flow  vs. 
the  along-shore  flow.  The  across-shore  flow  was  dominated  by  features  order  20  m  vertically  and 
about  2  days  in  length.  There  is  some  suggestion,  especially  at  ADCP  1 ,  that  these  features  were  propa¬ 
gating  upward.  Alongshore  currents  were  statistically  well  correlated  between  moorings  (r-.72-.86) 
but  across-shore  currents  were  weakly  correlated  (r**.3-.5).  The  origin  of  these  smaller-scale  across- 
shelf  motions  is  not  obvious  from  this  data  set.  The  smaller-scale  motions  may  have  an  alongshore 
component  as  well,  but  it  is  masked  by  the  larger-scale  alongshore  flows. 

4.  Model/data  comparisons 

4. 1.  Basic  description  of  the  three  numerical  models 

A  key  goal  of  the  ASAP  program  was  to  collect  data  sets  to  improve  numerical  models  for  pre¬ 
dicting  ocean  conditions  at  the  scales  necessary  to  resolve  coastal  upwelling  and  subsequent  plume 
dynamics.  The  models  were  run  in  near-real  time  to  provide  knowledge  input  to  the  adaptive  sampling 
decision  tree.  The  three  models  used  were  the  Harvard  Ocean  Prediction  System  (HOPS)  (Robinson, 
1999;  Haley  etal.,  2009;  Haley  and  Lermusiaux,  2010),  the  Jet  Propulsion  Laboratory  implementation 
of  the  Regional  Ocean  Modeling  System  (JPL/ROMS)  (Schepetkin  and  Williams,  2004;  Chao  etal.,  2009; 
Wang  et  al,  2009),  and  the  Navy  Coastal  Ocean  Model/Innovative  Coastal-Ocean  Observing  Network 
(NCOM/ICON)  model  (Shulman  et  al.,  2009,  2010).  Two  versions  of  the  HOPS  model  were  examined: 
HOPS  solo  (hereafter  HOPSs)  which  had  1.5  km  horizontal  resolution,  and  HOPS  nested  (hereafter 
HOPSn)  which  was  nested  within  HOPSs  and  had  0.5  km  horizontal  resolution  (Fig.  5).  JPL/ROMS  and 
NCOM/ICON  had  rectangular  grid  cells  near  the  moorings  with  about  1 .5  km  x  2.5  km  horizontal  res¬ 
olution.  All  the  models  had  30-32  sigma  coordinates  in  the  vertical.  While  conceptually  similar  in 
operation,  the  models  differed  significantly  in  terms  of  the  mixing  schemes  employed,  initial  condi¬ 
tions  used,  and  the  larger-scale  models  they  were  nested  within.  All  four  models  were  forced  using 
surface  wind  stresses  and  heat  fluxes  from  the  Coupled  Ocean/Atmosphere  Mesoscale  Prediction  Sys¬ 
tem  (COAMPS™)  (Doyle  et  al.,  2009),  and  assimilated  the  NPS  aircraft  SST,  glider  T  and  S,  dnd  AUVT 
and  S  data.  The  glider  data  set  was  one  of  the  largest  and  most  unique  data  sets  available  for  assim¬ 
ilation.  The  ten  gliders  adaptively  patrolled  the  22  x40kmx  1000  m  region  under  an  autonomous 
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Fig.  5.  Model  domains  off  central  California  for  the  four  numerical  models  compared  against  the  moored  ADCP  data  on  the 
continental  shelf.  The  mooring  positions  are  indicated  by  the  black  diamonds. 

coordinated  control  system  (Leonard  et  al.,  2010)  programmed  to  adhere  to  the  22  km  x  2.2  d  length 
and  time  correlation  scales,  respectively.  A  total  of  4530  Spray  profiles  and  10,619  Slocum  profiles 
were  collected.  These  data  did  a  good  job  of  resolving  the  eddies  and  cold  filaments  within  the  inten¬ 
sively  sampled  region.  Additionally,  the  JPL/ROMS  model  assimilated  T  and  S  data  from  the  MBAR1 
buoys  Ml  and  M2.  None  of  the  models  assimilated  velocity  from  the  vehicles,  moorings,  or  coastal  HF 
radar. 

The  models  all  use  different  mixing  parameterizations  (see  references  for  details).  The  effect  of  this 
on  mesoscale  model  output  differences  is  not  obvious,  but  is  limited  for  short  ocean  forecasts  of  a  few 
days  or  less.  Two  factors  that  clearly  are  important  are  the  model  initial  and  open  boundary  conditions. 
Given  the  small  scale  of  the  domain  off  Point  Ano  Nuevo,  the  boundary  forcing  is  extremely  important 
to  obtaining  the  correct  answer  within  the  model  domain.  Such  well-known  features  as  the  California 
Undercurrent  over  the  continental  slope  and  the  barotropic  poleward  flow  over  the  continental  shelf 
were  clearly  remotely  forced  and  entered  the  domain  through  the  open  boundaries.  The  JPL/ROMS 
model  had  the  longest  spin-up  time  using  climatology  since  July  2003,  using  boundary  conditions 
from  the  next-larger  ROMS  model  (5-km)  it  was  nested  within.  NCOM/1CON  likewise  used  boundary 
conditions  from  the  next  larger  nest,  called  NCOM  California  Current  System  (NCOM/CCS).  The  run 
was  initialized  using  more  recent  data  starting  July  11, 2006.  HOPS  solo  was  unique  among  the  three 
in  that  is  was  not  nested  into  a  larger-scale  model.  This  model  was  initialized  using  local  data,  that 
is,  the  early  glider  runs  and  ship  surveys  from  the  POINT  SUR  and  the  WECOMA  conducted  in  the 
area,  and  employed  modified  Orlanski  radiation  conditions  at  the  boundaries.  All  the  regional  models, 
either  directly  or  indirectly,  made  use  of  satellite  products  such  as  sea  surface  temperature  (SST)  and 
surface  height  (SSH).  Some  of  the  runs  (such  as  5  km  JPL/ROMS)  assimilated  satellite  data  directly, 
while  the  HOPS  runs  did  not  assimilate  SSH.  NCOM/CCS  pulled  in  SSH  via  an  intermediate  step,  the 
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Navy’s  Modular  Ocean  Data  Assimilation  (MODAS)  system.  This  information  was  thus  transferred  to 
the  smaller  1 .0-1 .5  km  nests  via  the  open  boundary  conditions. 

4.2.  Visual  comparisons 

To  develop  a  basis  for  comparison,  the  output  from  all  four  numerical  models  first  underwent 
the  same  low-pass  filtering  and  coordinate  rotation  as  the  ADCP  data.  The  principal  axis  ellipses  are 
interesting  in  their  own  right:  At  ADCP  1  the  NCOM  model  overestimated  the  currents  by  about  2  cm/s 
relative  to  both  the  remaining  models  and  the  data,  while  the  ROMS  model  output  was  misaligned  by 
12°  relative  to  the  other  estimators  (Fig.  6a).  At  ADCP  2,  all  the  models  statistically  underestimated 
the  true  currents  and  the  HOPSs  model  had  too  much  energy  in  the  across-shelf  component  (Fig.  6b). 

Three  depths  were  chosen  for  visual  and  statistical  comparison:  near-surface  ( 10, 1 2  m),  mid-depth 
(24,  52  m),  and  near-bottom  (46,  72  m)  for  ADCPs(1, 2),  respectively. 

Recall  that  at  ADCP  1,  the  observed  currents  consisted  of  strong  poleward  flow  from  3  to  8  August, 
followed  by  equatorward  flow  from  9  to  17  August,  and  weakly  poleward  flow  for  the  remainder 
of  the  record  (18-31  August).  This  was  true  at  all  depths  with  limited  vertical  shear  (Fig.  7).  The 
NCOM  model  fluctuations  appeared  to  lead  the  data,  most  noticeably  near  the  surface  (Fig.  7a).  After 
mid-August,  the  NCOM  model  overestimated  the  strength  of  the  poleward  flow  by  a  factor  of  2-4, 
worse  near  the  surface  than  near  the  bottom.  The  HOPSs  model  overestimated  the  poleward  flow 
from  1  to  3  August,  missed  the  equatorward  event  completely,  and  also  overestimated  the  second 
poleward  event  (Fig.  7a).  Like  the  NCOM  model,  the  model/data  agreement  improved  with  depth 
and  near  the  bottom  the  agreement  was  much  better  than  higher  in  the  water  column  (Fig.  7c).  The 
ROMS  model  showed  better  agreement  with  the  data  near  the  surface  than  the  NCOM  or  HOPSs 
models  but  still  had  difficulty  reproducing  the  equatorward  flow.  The  ROMS  currents  opposed  the 
observed  currents  during  9-13  August  and  were  too  weak  from  13  to  18  August.  The  ROMS  output 
was  the  least  sheared  among  the  models  and  looked  much  the  same  at  all  depths  (Fig.  7,  bottom 
panels).  The  HOPSn  model  showed  the  best  visual  agreement  with  the  data.  Like  HOPSs,  the  HOPSn 
time  series  began  with  too  much  poleward  flow,  perhaps  a  product  of  model  initialization  and  the 
short  spin-up  time  relative  to  the  other  two  models.  Once  it  got  on  track,  the  model  produced  the 
August  13-17  equatorward  event  and  had  close  to  the  right  magnitude  for  the  second  poleward 
event. 

The  observations  at  ADCP  2  were  more  vertically  sheared  than  at  ADCP  1.  At  ADCP  2,  the  near¬ 
surface  currents  were  entirely  equatorward  following  the  wind  stress  during  August  8-31 ,  while  the 
mid-depth  and  near  bottom  currents  were  poleward  from  August  18-28  in  opposition  to  the  surface 
flow  (Fig.  8,  compare  a  and  b).  The  observed  surface  currents  at  ADCP  2  also  opposed  the  surface 
currents  at  ADCP  1  from  August  18-28,  indicating  the  presence  of  a  velocity  front  between  the  two 
instruments  (compare  Figs.  7a  and  8a).  A  model’s  performance  will  therefore  be  quite  sensitive  to  its 
ability  to  correctly  position  that  front.  Like  at  ADCP  1,  the  NCOM  model  looked  good  for  the  first  half 
of  the  record,  again  leading  the  data  slightly,  but  deteriorated  during  the  second  half  of  the  month 
where  the  surface  poleward  flow  in  the  model  output  opposed  the  observed  flow  at  the  ADCP  ( Fig.  8a). 
This  was  perhaps  because  aircraft  and  glider  data  were  available  to  assimilate  during  August  1  -1 5  but 
not  during  August  16-31.  The  NCOM  model  and  data  agreement  at  ADCP  2  was  much  better  near 
the  bottom  where  the  model  captured  the  poleward/equatorward/poleward  transitions  although  the 
timing  was  still  off,  and  there  was  too  little  equatorward  flow  and  too  much  poleward  flow  at  the 
bottom.  The  equatorward  event  was  weak  in  the  ROMS  model  near  the  surface  and  absent  at  the  other 
depths,  indicating  it  misplaced  the  location  of  the  California  Current  boundary  (Fig.  8a-c).  It  also  badly 
overestimated  the  poleward  flow  after  August  1 8.  Like  the  NCOM  model,  the  currents  after  August  1 8 
more  closely  resembled  the  observations  at  ADCP  1  than  at  ADCP  2,  indicating  that  the  model  placed 
the  velocity  front  offshore  of  the  observed  velocity  front.  This  is  consistent  with  both  NCOM  and  ROMS 
currents  being  geographically  misplaced  a  bit  too  far  offshore  with  respect  to  the  observations.  The 
HOPSs  model  looked  good  at  the  beginning,  but  then  had  equatorward  rather  than  poleward  flow  at 
mid-depth  and  near-bottom  after  August  18  (Fig.  8b,  c).  The  HOPSn  model  looked  the  best  at  site  2  as 
it  did  at  site  1 .  It  faithfully  produced  both  the  opening  poleward  and  equatorward  flow  events,  and  at 
least  had  the  sign  correct  at  the  end. 
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Fig.  6.  Mid-water  column  principal  axis  ellipses  for  both  models  and  moorings  from;  (a)  ADCP  1 , 34  m  depth;  and  (b)  ADCP  2. 
54  m  depth.  North  is  towards  the  top  of  the  plot. 
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Fig.  7.  A  comparison  of  currents  from  ADCP  1  and  the  four  numerical  models  for:  (a)  near-surface  ( 1 0  m),  (b)  mid-depth  (24  m). 
and  (c)  near-bottom  (46  m).  Poleward  alongshore  is  up. 
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Fig.  8.  A  comparison  of  currents  from  ADCP2  and  the  four  numerical  models  for:  (a)  near-surface  (12  m),(b)  mid-depth  (52  m), 
and  (c)  near-bottam  (72  m). 


43.  Statistical  comparisons 

While  already  interesting,  these  stick  vector  plots  form  merely  an  introduction  to  more  quan¬ 
titative  statistical  comparisons.  The  tools  used  to  facilitate  these  comparisons  were  the  temporal 
cross-correlation  functions  and  the  root  mean  square  (RMS)  error  between  the  model  output  and 
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the  ADCP  data  at  several  grid  points  surrounding  each  ADCP  location  (Fig.  9).  The  correlations  deter¬ 
mine  if  the  model  output  was  geographically  or  temporally  displaced  from  the  data,  while  the  RMS 
error  estimates  the  magnitude  of  the  model/data  differences. 

For  the  closest  model  grid  point  to  the  ADCP  location  and  each  of  the  surrounding  grid  points,  the 
correlation  coefficient  pxy( r)  was  calculated  as  (Bendat  and  Piersol.  1986): 


where  Cxy(r)  is  the  cross-covariance  between  the  model  output  (x)  and  the  data  (y),  oXt  ay  are  the 
variances  of  series  x  andy,  and  r  is  the  time  lag.  Significance  levels  for  the  correlations  were  calculated 
following  the  method  of  Amos  and  Koopmans  (1963)  using  the  effective  degrees  of  freedom  estimates, 
Neff,  defined  as  (Davis,  1976): 


where  N  is  the  length  of  the  time  series  in  hours  and 


00 


Tn=  ^C0(iAt)Cn(iAt)At 


is  the  correlation  time  scale  for  the  data.  More  self-correlated  data  results  in  fewer  true  degrees  of 
freedom  and  therefore  higher  significance  levels  for  the  correlation  coefficient. 

When  examining  the  results  from  these  calculations  it  is  important  to  keep  the  spatial  distribution 
of  the  points  being  compared  in  mind  since  they  were  different  for  each  model  (Fig.  9).  For  the  NCOM 
and  both  HOPS  models,  the  grid  points  were  oriented  along/across-shore  while  the  ROMS  model  grid 
was  oriented  north/south.  The  HOPSn  grid  was  about  3  times  more  resolved  than  the  others  and 
less  variation  is  thus  expected  in  the  correlation  coefficients  between  these  points  and  the  ADCP  at 
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Fig.  9.  Grid  point  locations  used  to  facilitate  the  spatial  correlations,  with  the  ADCP  located  at  the  center  of  each  grid  as  indicated 
by  the  asterisk. 
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the  center.  The  results  (Appendix  A)  are  displayed  as  box  charts  showing  the  correlation.  RMS  error, 
and  optimum  lag  time  for  each  of  the  nine  grid  points  compared  to  the  ADCP  located  at  the  center. 
Statistically  insignificant  results  are  shaded,  and  time  lags  are  presented  only  for  the  statistically 
significant  results.  A  negative  (positive)  phase  means  the  model  leads  (lags)  the  data.  An  empty  box 
indicates  no  data,  usually  because  the  selected  model  depth  was  deeper  than  the  bottom  on  the  inshore 
side  of  the  mooring.  For  display,  the  along-  and  across-shore  components  for  each  model  from  the  grid 
point  closest  to  the  ADCP  were  compared  with  the  data  ( Fig.  1 0).  Then,  the  tables  were  used  to  produce 
component  plots  of  the  “best"  model/data  comparisons  for  each  model  using  the  grid  point  with  the 
highest  correlation  shifted  relative  to  the  data  by  the  corresponding  time  lag  (Fig.  1 1).  A  few  robust 
results  emerged  from  these  comparisons. 

For  the  alongshore  component,  all  the  models  except  NCOM  at  ADCP  1  were  significantly  correlated 
with  the  data.  At  site  1,  the  ROMS  model  produced  the  highest  correlations  followed  by  HOPSs  and 
HOPSn.  The  ROMS  model  output  was  very  nearly  in  phase  with  the  observations  while  both  the  HOPS 
models  lead  the  data  by  about  1 6-24  h  (Fig.  10,  left  column).  At  site  2,  the  highest  correlations  were 
at  HOPSn  followed  closely  by  HOPSs  and  ROMS  where  the  correlation  levels  were  quite  similar.  The 
HOPSn  model  lead  the  data  by  about  5  h  at  the  surface  but  lagged  by  16-21  h  subsurface  The  HOPSs 
model  also  lead  the  data  but  by  smaller  amounts,  0-8  h  for  the  most  correlated  points.  The  phase  num¬ 
bers  for  ROMS  at  site  1  werea  bit  erratic:  Choosingthe  most  correlated  points,  the  model  was  in  phase  at 
the  surface,  lead  the  data  by  22  h  at  mid-depth,  and  lagged  by  7  h  at  the  bottom  (Appendix  Al.  In  terms  of 
the  geographic  accuracy,  HOPSn  at  site  2  was  consistently  better  correlated  at  the  NE  point  rather  than 
at  the  center  point,  perhaps  indicating  a  slight  misalignment  of  model  features,  although  the  distances 
were  small.  None  of  the  other  models  did  consistently  better  or  worse  at  the  other  model  grid  points. 

The  HOPS  nested  model  was  the  only  one  with  any  skill  at  reproducing  the  across-shelf  component. 
Some  significant  correlations  were  observed  for  the  ROMS  model  but  they  appeared  fortuitous  since 
they  were  180  out  of  phase.  At  site  2,  the  correlations  between  HOPSn  and  the  data  were  well  above 
the  significance  levels  (Appendix  A)  and  the  visual  agreement  with  the  observations  was  good  (not 
shown).  This  indicates  that  HOPSn  responded  well  to  the  local  wind  forcing  as  described  earlier  from 
the  observations  alone,  which  was  more  evident  at  site  2  than  at  site  1.  The  HOPSn  across-shore 
component  always  lagged  behind  the  observations  by  longer  (18  h)  at  mid-depth  than  al  the  surface 
and  bottom  (8  h).  This  is  again  consistent  with  a  boundary  layer  response,  however  it  takes  too  long  in 
the  model  for  the  momentum  to  be  transferred  to  mid-depth.  It  seems  that  very  high  resolution  models 
are  necessary  to  successfully  simulate  the  across-shore  component  of  the  flow  on  the  continental  shelf 
in  this  region.  This  is  likely  due  to  the  shorter  correlation  scales  in  the  across-shelf  direction  than  the 
alongshelf  direction  (Ramp  et  al.,  2008). 

5.  Discussion 

5.  /.  Remote  forcing,  data  assimilation,  and  boundary  conditions 

All  three  models,  with  well-established  performance  at  larger  space  and  time  scales,  had  difficulty 
reproducing  the  current  variability  in  this  small  sample  region  off  Point  Ano  Nuevo.  One  potential 
explanation  is  that  these  small  domains  cannot  correctly  capture  remote  forcing  in  the  form  of  along¬ 
shore  pressure  gradient  forces  or  coastally  trapped  waves,  which  propagate  from  south  to  north  with 
the  coast  on  the  right  in  this  region.  This  means  that  larger-scale  features  such  as  these  must  be 
captured  by  the  assimilation  scheme  or  fed  to  the  smaller  regional-scale  models  via  their  boundary 
conditions. 

Considering  data  assimilation  first,  the  circulation  models  in  ASAP  were  quantitatively  compared 
with  the  observed  currents  at  the  two  ADCP  moorings  located  to  the  south  of  the  area  that  was  inten¬ 
sively  sampled  by  the  gliders  and  aircraft.  Under  this  configuration,  poleward-propagating  events  of 
any  kind  reached  the  ADCPs  ahead  of  the  major  sampling  area  around  the  upwelling  center.  The  lack 
of  data  for  assimilation  to  the  south  of  the  ADCPs  is  one  of  the  reasons  why  the  models  represented 
the  effects  of  the  remote  forcing  differently.  All  of  the  real-time  models  used  sequential  data  assimi¬ 
lation  schemes  which  do  not  correct  initial  and  boundary  fields  backward  in  time.  In  such  a  situation, 
the  boundary  and  initial  conditions  chosen  by  the  modelers,  the  data  and  background  covariances 
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Fig.  10.  Model/data  (gray/black  lines)  comparisons  for  the  alongshore  current  components  at  ADCP  1  (left  column)  and  ADCP 
2  (right  column).  The  three  vertical  panels  for  both  columns  are  the  comparisons  for  (a)  near-surface  (10,  1 2  m),  (b)  mid-depth 
(24.  52  m).  and  (c)  near-bottom  (46.  72  m). 
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Fig.  11.  As  in  Fig  1 0.  except  the  plots  are  for  the  grid  points  (see  Fig.  9)  and  time  lags  that  resulted  in  the  best  agreement  between 
the  model  output  and  the  ADCP  data,  as  judged  by  the  highest  correlation  coefficient  (Appendix  A).  The  notation  indicating  the 
position  of  the  winning  grid  point  relative  to  the  ADCP  is  the  same  as  points  on  a  compass  with  a  "C  indicating  the  center  point 
(no  geographical  offset).  The  convention  for  the  phase  is  that  a  negative  (positive)  sign  means  the  model  leads  (lags)  the  data. 
For  example  (upper  left  panel),  the  alongshore  component  of  the  ROMS  model  at  10  m  showed  the  best  agreement  with  the 
data  when  the  southern  grid  point  was  lagged  6  h  in  time.  The  gray  shading  indicates  panels  where  the  correlation  coefficient 
was  not  statistically  significant,  but  the  phase  shift  was  used  anyway  for  completeness. 
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Fig.  1 2.  Sea  surface  height  (SSH)  anomalies  from  24  N  to  40  N  as  computed  by  (left  panel)  the  Navy  regional  NCOM  CCS  model 
with  global  NCOM  boundary  conditions  and  (right  panel)  the  global  HYCOM  model.  The  observed  SSH  from  coastal  sea  level 
observations  at  Monterey  (36  36'N)  and  San  Diego  (32  43'N)  are  included  as  the  heavy  black  lines.  Time  series  at  the  same 
locations  sub-sampled  from  the  two  different  model  configurations  are  shown  as  the  red  (NCOM)  and  blue  (HYCOM)  lines, 
respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the 
article.) 


used  in  the  assimilation,  and  the  nesting  schemes  as  well  as  the  model  resolutions  contributed  to  the 
differences. 

Two  of  the  models,  ROMS  and  NCOM,  obtained  their  boundary  conditions  from  the  models  they 
were  nested  within.  The  boundaries  in  HOPSn  were  obtained  from  HOPSs,  which  was  initialized  using 
the  ASAP  opening  ship  surveys  and  other  available  data.  The  HOPSs  boundary  conditions  were  subse¬ 
quently  updated  using  fits  of  the  model  to  all  available  data.  To  explore  the  effects  of  different  boundary 
conditions  and  coordinate  systems,  additional  model  runs  were  made  subsequent  to  the  ASAP  pro¬ 
gram  using  the  Navy’s  global  Hybrid  Coordinate  Ocean  Model  (HYCOM)  to  force  the  NCOM/1CON 
model,  rather  than  the  global  NCOM  via  the  NCOM  CCS  output.  Forcing  the  same  NCOM  ICON  model 
directly  from  the  global  HYCOM  model  produced  a  dramatic  improvement  in  the  model/data  compar¬ 
isons:  The  complex  correlations/angular  displacements  (Kundu,  1976)  improved  from  0.04/52.0°  to 
0.69/- 1 1 .5°  at  ADCP 1 .  and  from  0.1 7/2.9°  to  0.86/  *4.9°  at  ADCP  2.  This  difference  represents  a  change 
from  essentially  no  predictive  skill  to  very  good  predictive  skill. 

These  improvements  can  be  understood  by  comparing  the  sea  surface  height  (SSH)  anomalies 
along  the  California  coast  calculated  from  the  NCOM/NCOM  CCS  model  and  the  HYCOM  n)odel  with 
the  coastal  SSH  observations  (Fig.  12).  Both  models  and  the  observations  show  that  the  SSH  anomalies 
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during  August  2006  were  larger  than  normal  and  much  larger  than  during  AOSN-11  in  summer  2003 
(Ramp  et  al.,  2009;  Shulman  et  al.(  2010).  Both  models  also  show  that  these  SSH  anomalies  were  propa¬ 
gating  poleward  at  1  -3  m  s  1 ,  as  can  be  seen  by  the  slope  of  the  maxima  (red  region)  in  the  latitude  vs. 
time  plots  (Fig.  12).  These  anomalies  therefore  represent  remote  forcing  which  will  impact  the  smaller 
ASAP  region  off  Point  Afio  Nuevo.  Comparing  the  model  SSH  (red,  blue  lines)  with  the  observed  SSH 
at  Monterey  and  San  Diego  (black  lines)  shows  that  NCOM/NCOM  CCS  model  underestimated  these 
propagating  anomalies  with  respect  to  both  the  observations  and  the  HYCOM  model,  which  closely 
tracked  the  observed  fluctuations. 

Many  previous  authors  have  documented  poleward-propagating  motions  along  the  U.S.  west  coast 
in  coastal  sea  level  records  (e.g.  Chelton  and  Davis,  1982;  Denbo  and  Allen,  1987;  Spillane  et  al.,  1987) 
current  meter  records  (Ramp  et  al.,  1997a)  and  model  output  (Pares-Sierra  and  O’Brien,  1989;  Ramp 
et  al.,  1997b).  For  ASAP  2006,  a  coherence  and  phase  analysis  of  the  observed  SSH  anomalies  relative 
to  the  annual  mean  along  the  coast  during  May  to  October  2006  showed  poleward  propagation  at 
2.7  m  s  1  between  Port  San  Luis  and  Monterey  for  oscillations  with  a  15-day  period.  Using  a  simple 
c  =  w/fc  relationship  where  c  is  the  phase  speed,  co  is  the  angular  frequency,  and  k  is  the  horizontal 
wave  number,  these  fluctuations  would  have  a  wavelength  of  3695  km  if  they  were  propagating  as 
free  waves  along  the  coast.  Previous  authors  found  waves  along  the  Big  Sur  coast  in  the  29-day  band 
that  closely  resembled  linear,  coastally  trapped  mode-1  internal  Kelvin  waves  (Ramp  et  al.,  1997a). 
The  dispersion  relation  for  these  waves  is  given  by 

2  N2k2 

w  = - r- 

k2  +  (rm/H)2 

where  N  is  the  Brunt-Vaisala  frequency,  n  is  the  mode  number,  H  is  the  bottom  depth,  and  the 
other  variables  are  as  previously  defined.  Choosing  N  =  5.5  x  10  3  s  1  as  a  representative  value  and 
H=  1500  m  as  a  mean  depth  over  the  slope,  the  dispersion  relation  gives  a  wavelength  of  3490  km 
for  the  15-day  period.  The  fluctuations  observed  in  both  coastal  sea  level  and  the  Navy  models  are 
therefore  consistent  with  coastally  trapped  mode-1  internal  Kelvin  waves.  It  appears  likely  therefore 
that  the  improvement  in  the  model/data  comparisons  at  ADCP  1  and  ADCP  2  with  HYCOM  nesting 
were  due  to  that  model’s  more  accurate  representation  of  the  poleward  fluctuations  entering  the 
ASAP  region  from  the  south.  Additional  model  runs  (not  shown)  changing  only  smoothed  (sigma)  vs. 
un-smoothed  (hybrid)  bottom  coordinates  showed  improvements  using  the  hybrid  coordinates,  but 
these  improvements  were  much  smaller  than  those  resulting  from  using  the  HYCOM  vs.  NCOM  model. 
This  indicates  that  the  improved  boundary  conditions  were  the  most  important  factor  in  improving 
the  NCOM  ICON  model  performance. 


5.2.  Importance  of  offshore  eddies 

Given  the  relatively  good  visual  and  statistical  comparisons,  a  few  plan-view  plots  from  the  HOPSn 
model  were  chosen  to  illuminate  the  larger-scale  context  for  the  ADCP  observed  currents  (Fig.  13). 
The  first  pair  on  August  4  and  6  (Fig.  13a  and  b)  show  the  development  of  the  strong  poleward  flow, 
one  of  the  most  prominent  features  of  the  record.  The  cold  water  at  Point  Ano  Nuevo  was  quickly 
swamped  by  the  warm  water  from  the  south  during  this  event.  The  only  really  strong  equatorward 
flow  during  August  2006  took  place  during  the  days  surrounding  August  11.  This  was  due  to  the 
development  of  a  coastal  upwelling  jet  between  the  upwellingcenterand  a  warm  shoreward  incursion 
of  the  California  Current  offshore  (Fig.  13c).  During  subsequent  upwelling  events,  as  exemplified  for 
instance  by  conditions  on  August  23  (Fig.  13d),  the  upwelling  front  was  located  further  offshore  and 
the  equatorward  flow  at  the  two  ADCPs  was  weaker  or  non-existent.  A  comparison  of  Fig.  13a  and  c 
supports  the  conclusion  of  earlier  investigators  that  the  bifurcation  of  the  upwelling  plume  is  closely 
related  to  the  position  and  strength  of  the  offshore  eddy  field.  When  the  eddy  was  well  offshore 
(Fig.  13a)  the  cold  plume  likewise  advected  offshore.  When  the  eddy  was  strong  and  closer  to  shore 
(Fig.  13c)  the  cold  plume  moved  primarily  southward. 
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Fig.  1 3.  HOPSn  model  currents  ( black  vectors)  overlaid  on  temperature  (color  bar)  for  the  sea  surface  layer  off  Point  Aflo  Nuevo 
for  (a)  4  August,  (b)  6  August,  (c)  11  August,  and  (d)  23  August  2006.  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


6.  Summary  and  conclusions 

The  ocean’s  response  to  both  upwelling  and  relaxation  events  was  studied  for  the  month  of 
August  2006  near  Point  Ado  Nuevo,  CA  as  part  of  the  ASAP  MB06  field  experiment.  An  aircraft, 
several  autonomous  gliders  and  AUVS,  moored  ADCPs,  and  three  numerical  models  were  used  to 
describe  the  flow.  This  paper  provides  a  basic  description  of  the  flow  from  the  aircraft  and  moor¬ 
ings  and  compares  those  results  with  numerical  simulations  from  the  HOPS,  JPL/ROMS,  and  NCOM 
models. 

The  aircraft  produced  a  high-resolution  spatial  description  of  the  wind  stress  variability  which 
was  previously  unavailable  for  the  region.  Several  new  features  emerged  including  the  wind  stress 
maximum  to  the  north  of  the  Point  during  poleward  wind  events,  a  clear  maximum  offshore  from  the 
Point  during  upwelling  events,  and  a  strong  wind  stress  curl,  also  during  upwelling  favorable  events, 
which  could  potentially  drive  upwelling  with  strength  (5-10  md  *)  equal  to  the  coastal  divergence. 
This  spatial  variability  in  the  wind  itself  likely  plays  an  important  role  in  determining  the  location 
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of  3-D  upwelling  centers  along  the  west  coast  of  the  United  States.  The  wind  stress  observed  at  the 
closest  local  buoy  was  often  not  representative  of  the  actual  wind  stress  off  the  Point  or  to  the  north 
of  the  Point  and  much  caution  is  needed  when  applying  the  buoy  wind  observations  to  the  local 
coastal  oceanography.  The  MBARI  buoy  M2  was  nevertheless  the  most  representative  time  series 
available  for  the  entire  month  of  August  since  the  aircraft  was  not  available  after  August  15.  The 
buoy  recorded  upwelling-favorable  wind  events  with  approximately  equal  strength  during  August 
1-3,  8-12,  16-25,  and  28-31.  Poleward  winds  less  than  2.5 ms-1  blew  during  August  4-7,  13-15, 
and  26-27. 

Despite  being  only  6.5  km  apart,  moorings  1  and  2  responded  differently  to  the  sequence  of 
upwelling  favorable  winds  separated  by  brief  relaxations.  The  flow  at  site  1  was  almost  always  pole- 
ward  and  often  appeared  to  surge  out  of  the  Monterey  Bay.  Site  2  was  just  far  enough  offshore  to  miss 
this  phenomenon.  The  site  2  flow  was  more  strongly  sheared  vertically  and  more  coherent  with  the 
local  wind  stress,  especially  the  across-shore  current  component  with  the  alongshore  wind.  Still,  there 
were  many  times  when  the  currents  at  site  2  were  not  coherent  with  the  winds  and  were  instead  well 
correlated  with  the  position  of  the  mesoscale  features  in  the  coastal  ocean.  The  strongest  equatorward 
currents  were  clearly  associated  with  the  upwelling  jet,  and  the  across-shore  position  of  the  jet  largely 
determined  the  current  strength  at  site  2,  less  so  at  site  1  which  was  located  well  inshore  of  the  jet 
most  of  the  time. 

Given  the  very  small  scales  of  motion  and  the  strong  shears  between  the  two  moorings,  it  is 
perhaps  not  surprising  that  the  three  models  all  had  difficulty  simulating  the  flow.  The  models 
had  more  success  during  the  first  half  of  the  month  when  the  flow  was  more  wind-driven  and 
there  were  aircraft  data  available  to  assimilate.  The  models  struggled  during  the  second  half  of 
the  experiment  when  the  flow  appeared  to  be  more  remotely  forced  and  local  observations  were 
lacking.  A  post-experiment  analysis  of  observed  sea  level  anomalies  along  the  coast  and  output 
from  the  global  HYCOM  model  re-analysis  both  showed  that  sea  level  fluctuations  during  August 
2006  were  anomalously  strong  and  propagated  poleward  with  characteristics  consistent  with  free 
coastally  trapped  mode-1  internal  Kelvin  waves.  The  failure  of  the  outer-nests  driving  the  open 
boundaries  of  the  JPL-ROMS  and  NCOM  ICON  models  to  correctly  capture  these  fluctuations  may 
account  for  many  of  the  observed  weaknesses  in  the  model/data  comparisons.  Subsequent  re- 
analyses  of  the  NCOM  ICON  model  forced  by  the  global  HYCOM  model,  which  more  accurately 
represented  the  coastal  sea  level  variability,  showed  a  dramatic  improvement  in  model  predictive 
skill. 

The  only  model  that  had  any  demonstrated  skill  simulating  the  across-shore  component  was  the 
HOPS  nested  model  (HOPSn),  which  had  a  three  times  finer  grid  resolution  than  the  others  (0.5  vs. 
1.5  km).  This  very  fine  simulation  scale  appears  to  be  necessary  to  correctly  place  the  strong  fronts 
that  often  occur  in  the  coastal  ocean. 
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Appendix  A. 

Statistical  comparisons  of  the  model  output  at  nine  grid  points  with  the  ADCP  time  series,  co¬ 
located  at  the  center  of  the  box.  The  maximum  correlation  coefficient  at  the  optimum  lag  (h),  and 
RMS  error  between  instruments  are  shown  for  each  comparison.  Correlations  that  are  not  statistically 
significant  have  been  grayed  out.  The  across-  and  alongshore  components  are  shown  in  the  left  and 
right  column,  respectively.  The  surface,  mid-depth,  and  near-bottom  time  series  are  at  the  top,  middle, 
and  bottom. 
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